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Abstract 



The stochastic inflation program is a framework for understanding the dy- 
namics of a quantum scalar field driving an inflationary phase. Though widely 
used and accepted, there have over recent years been serious criticisms of this 
theory. In this paper I will present a new theory of stochastic inflation which 
avoids the problems of the conventional approach. Specifically, the theory 
can address the quantum-to-classical transition problem, and it will be shown 
to lead to a dramatic easing of the fine tuning constraints that have plagued 
inflation theories. 
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I. INTRODUCTION 



The inflationary universe scenario asserts that, at some very early time, the universe 
went through a de Sitter phase expansion with scale factor a(t) growing as e Ht . Inflation 
is needed because it solves the horizon, flatness and monopole problems of the very early 
universe and also provides a mechanism for the creation of primordial density fluctuations. 
For these reasons it is an integral part of the standard cosmological model |IJ . 

The inflationary phase is driven by a quantum scalar field with a potential V(<E>), that 
can take on many different forms that satisfy the 'slow roll' conditions. Inflationary scenarios 
fall into two types. In the first, which includes the so-called old and new inflation, the scalar 
field driving inflation (the inflaton) is assumed to be in thermal equilibrium with the rest of 
the universe. The universe is assumed to obey the standard hot big bang cosmology in the 
period preceding and after inflation. In the chaotic inflation scenario these assumptions are 
relaxed. Instead the inflaton is assumed to be only very weakly coupled to other fields. This 
makes it possible to choose initial states for the inflaton which are far from equilibrium. The 
standard big bang cosmology now only applies after the reheating stage of inflation. 

In the conventional approach to inflaton dynamics [Q, the inflaton field $ is first split 
into a homogeneous piece and an inhomogeneous piece 

$(x,s)=0(s)+^(x,s), (1.1) 

where <fi is interpreted as the field $ coarse-grained over a volume Q: 

#0 = ^ / $(x,s)d 3 x. (1.2) 

In this description all information on scales smaller than the coarse-graining volume, such 
as density fluctuations, will be contained in if>. This means the field if) has a low frequency 
cutoff in its spectrum. The coarse-graining volume is set by the causal horizon which leads 
to Q ~ H~ 3 . We will refer to this (f> as the global order parameter since it is constructed from 
a coarse-graining volume that is always larger than the observable universe. The dynamics 
of the global order parameter is then postulated to obey the classical 'slow roll' equation of 
motion 

^ = 0. (1.3) 

This equation governs the dynamics of (f> which drives the inflationary phase. It is also 
possible to discuss the generation of primordial density fluctuations using if). Assuming that 
<f) ^> if), it can be shown that if> is described by a free massless minimally coupled quantum 
scalar field. During exponential inflation the quantum fluctuations of if) grow as (if) 2 ) ~ H 3 t. 
These quantum fluctuations are then identified with the classical fluctuations which generate 
primordial density fluctuations [0,0. 

Although the conventional theory sketched above is widely accepted, it must be consid- 
ered problematic for three important reasons: 

• it is assumed, without justification, that the global order parameter <f> can be treated 
as a classical order parameter, and that the quantum fluctuations of if) are equivalent 
to classical fluctuations. 
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• the fluctuations of ip lead to an overproduction of primordial density fluctuations. This 
can only be avoided by unnaturally fine tuning the coupling constants in the inflaton 
potential. 

• the backreaction of ip on the dynamics of is ignored. 

Based on the conventional theory there should be no reason to expect <fi and ip to behave 
classically. The action for is non-linear and strongly time dependent. Both of these 
properties will always generate quantum mechanical coherence over time. Similarly, because 
ip is treated as a free quantum field its quantum vacuum state will always remain spatially 
homogeneous. It is only when these fluctuations become classical that spatial inhomogeneity 
is generated. Since the conventional theory treats and ip as independent closed systems 
it is impossible in principle for it to explain the quantum-to-classical transition of and ip. 
The purpose of this paper is to develop a more considered approach to inflaton dynamics 
which can address these three critical issues. 

The problems outlined above can all be addressed by a careful application of the princi- 
ples and techniques of non-equilibrium quantum statistical physics. When we coarse-grain 
the scalar field as in equation (1.2), what we are effectively doing is splitting our closed 
quantum system $ into a system sector and an environment sector ip. The system sector 
includes only wavelengths long compared to the coarse-graining scale while the environment 
sector comprises wavelengths short compared to the coarse-graining scale. By discarding 
information about the environment sector ip, and assuming that the system and environ- 
ment interact, we turn the unitary, reversible dynamics of the closed system $, into the 
nonunitary and irreversible dynamics of an open system 0. Thus we have gone from a fun- 
damental theory to an effective theory. The backreaction of any environment sector on a 
system sector will generally manifest itself by introducing noise, dissipation and renormal- 
ization effects into the dynamics of the system || . The noise and dissipation will be related 
by some fluctuation-dissipation relation. Noise will induce diffusion and decoherence (the 
loss of quantum mechanical coherence). Decoherence is the critical ingredient if we are to 
dynamically demonstrate the quantum-to-classical transition of an open system. 

In this approach our only interest in the environment ip is its backreaction on the system 
dynamics. In this case the global order parameter is no longer a suitable system variable 
because it contains no information about spatial structure within our universe. A system 
variable that does contain this information can be obtained simply by choosing the apparent 
rather than the causal horizon as the coarse-graining scale. This corresponds to a coarse- 
graining volume of Q ~ (He Ht )~ 3 . With this averaging volume we will refer to as the local 
order parameter. The use of a local rather than global order parameter for inflaton dynamics 
has been strongly advocated by Morikawa @]. In this scheme it is the classical fluctuations 
5(j) of the local order parameter which lead to density fluctuations rather than the quantum 
fluctuations derived from (ip 2 ). The new role of the field ip is to provide a noise source 
(via backreaction) in the quantum dynamics of the local order parameter. This noise will 
generate quantum decoherence which is the process that creates the classical fluctuations 
5(p of the local order parameter. This fundamental conceptual shift regarding the role of 
ip is the key to developing a theory that can address the quantum-to-classical transition 
problem. As a bonus it will be shown that the new theory leads to a dramatic easing of the 
fine tuning constraints that plague the conventional approach to inflaton dynamics. 
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A first attempt to study the backreaction of ip on the local order parameter 0, can be 
found in the 'stochastic inflation' program initiated by Starobinsky H and further developed 
by others ||. This program claims that the semiclassical equation of motion for the local 
order parameter is given by the Starobinsky equation 

• V'(<b) H 3 / 2 

where F w (t) is a zero mean gaussian white noise source of unit amplitude. The system 
field equation is thus transformed into a classical Langevin equation with a white noise 
source. A Fokker-Planck equation can also be derived which depicts the evolution of the 
probability distribution of the scalar field. Since it was first derived the Starobinsky equation 
has stimulated numerous studies. Much effort has been devoted to the solution of this 
stochastic equation for a description of the inflationary transition and the generation of 
primordial density fluctuations [0]. The same equation forms the theoretical foundation for 
studies on the very large scale structure of the universe ||. It has been claimed that the 
stochastic inflation program can explain the quantum-to-classical transition of the coarse- 
grained field |J. A quantum mechanical description, for which the Starobinsky equation is 
the semiclassical limit, has been developed using methods to quantize classically dissipative 
systems [|1(| (this list of references is by no means complete). 

In the derivation of the Starobinsky equation, interactions between the local order pa- 
rameter and its short wavelength environment i/j are derived from the quadratic terms in 
the scalar field action. Interactions that derive from a self-interacting potential are assumed 
to be small and are neglected. This procedure cannot be justified because, for free fields, 
there is no mode-mode coupling and therefore no way for short wavelengths to backreact 
on long wavelengths. Mode-mode coupling can only be generated by a self-interacting po- 
tential. Advocates of stochastic inflation would say that the coupling is generated by the 
time dependent nature of the system environment split. However it is simply impossible 
for a physical coupling to be generated by a time dependent choice of what constitutes the 
system and environment [ft^, 16|, lq| . We would therefore expect there to be, as long as self- 



interaction is ignored, no coupling and hence no backreaction of the environment ip on the 
local order parameter 0. This implies that classical fluctuations in the local order parameter 
must vanish for a free theory. We therefore conclude that the Starobinsky equation fails 
to describe the backreaction of ip on the local order parameter 0. 

A better interpretation for the Starobinsky equation is that it models the conventional 
approach to inflaton dynamics with a classical stochastic dynamical system. The left hand 
side is the usual postulated slow roll equation of motion for the global order parameter. The 
stochastic term is meant to classically model the quantum fluctuations of ip. This interpre- 
tation is plausible because solving the free Starobinsky equation gives (0 2 ) ~ H 3 t. This is 
the same as the the growth of quantum fluctuations in ip in the conventional approach. How- 
ever, since it is essentially equivalent to the conventional approach, Starobinsky's stochastic 
inflation program suffers from the same serious problems that were previously outlined. 

Several authors have previously made similar criticisms of Starobinsky's stochastic infla- 
tion and the conventional inflaton dynamics. Hu and Zhang |12| first questioned the ability 
of a time- dependent split in generating noise for a free field, and introduced the coarse- 
grained effective action for a proper treatment of backreaction. Habib (IT3] demonstrated 
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that the Starobinsky theory of stochastic inflation does not address the quantum-to-classical 
transition despite claims that it does |J. Hu, Paz and Zhang |16| (see also |17|JT8| ) provided 



an extensive critique of the Starobinsky theory and developed their own theory of stochastic 
inflation for the case of a quartic potential. Morikawa |l3j pointed out that self-interaction 
is necessary for the generation of fluctuations, and that a lack of understanding of this point 
is the cause of the fine tuning problem of inflation. This point has also been emphasized 
by Calzetta and Hu [I5| in a recent discussion on the problems with conventional inflaton 
dynamics. See also Morikawa |4| for a general discussion on inflaton dynamics. 

Clearly a new theory of stochastic inflation is required which can address the 3 problems 
of the conventional approach to inflaton dynamics. This new theory will require a fully 
nonequilibrium quantum mechanical formalism which deals with the statistical nature of 
mixed states and the dynamics of reduced density matrices. An ideal formalism for this 
is the Feynman Vernon influence functional The influence functional describes the 



averaged effect of the environmental degrees of freedom on the system degrees of freedom 
to which they are coupled. With the influence functional we can unambiguously identify 
a noise and dissipation kernel related by some generalized fluctuation-dissipation relation. 
The formalism leads directly to a propagator for the reduced density matrix which can then 
be used to derive a master equation. The propagator can be used to study the decoherence 
process which is the critical part of the quantum-to-classical transition. The master equation 
can be transformed into a Fokker-Planck type equation for the Wigner function. In the 
semiclassical limit we have a generalized Langevin equation. 

The stochastic inflation program, and indeed many other problems in the early universe 
Hj| , fall under the paradigm of quantum Brownian motion (QBM). QBM is one of the two 



major paradigms of non-equilibrium quantum statistical mechanics which is amenable to 
detailed analysis (the other being Boltzmann's kinetic theory). The complexity of the prob- 
lems in cosmology led Hu, Paz and Zhang p0| to first consider a systematic study of QBM in 
a general environment using the influence functional formalism. The special features associ- 
ated with a nonohmic bath, or ohmic bath at low temperatures are the appearance of colored 
noise and nonlocal dissipation. The methodology and viewpoint of QBM have been applied 
to the analysis of some basic issues in quantum cosmology |25,I8]21-27[|, effective field theory 
28| , |29[| , and the foundation of quantum mechanics, such as the uncertainty principle |]30|j3T|| 



and, most significantly, decoherence [p2|-p6| in the quantum-to-classical transition problem. 
(See the reviews of |37H39| and references therein for the standard literature on this topic). 
The QBM models in [^U] were further generalized by Hu and Matacz j|D| by making the 
system and bath oscillators the most general time dependent quadratic oscillators, and by 
considering more general system-bath couplings. The time dependence generates parametric 
amplification (squeezing) in the oscillators which is precisely how quantum fluctuations of a 
free field are amplified in the inflationary universe ||42|| . This formalism allows for a general 
study of nonequilibrium quantum statistical processes in time dependent backgrounds, as is 
the case in the early universe. 

Using the same methodology as in the previous QBM studies, Hu, Paz and Zhang \ W{ 



were the first, based on the system environment split orginally proposed by Starobinsky, 
to correctly formulate the problem of the backreaction of short wavelength quantum fluc- 
tuations ip, on the local order parameter 0. For more recent work of a similar nature see 
Lombardo and Mazzitelli These works are based on a rigorous quantum field theory of 
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open systems, however they do have some practical shortcomings. The dynamical systems 
obtained are of a functional nature further complicated by non-local dissipation and colored 
noise. It is very difficult to analyse the quantum and semiclassical dynamics of these systems. 
The models are restricted to flat space or conformally coupled fields in a de Sitter phase. 
For applications to inflation we are mainly interested in minimally coupled scalar fields in a 
de Sitter phase. The models are very specific, involving a A<3? 4 self-interaction. Deriving the 
renormalized effective action for the long wavelength sector is a complex calculation that 
can only be done in the context of specific potentials. In inflation there is great interest 
in a wide variety of potentials. In this case the quantum field theory calculations may be 
much more difficult to implement than for the A<3? 4 case. Because of these problems there 
has not yet been any detailed investigation of inflaton dynamics based on these quantum 
field theoretic open systems. 

In order to proceed further we need to develop a new simplified theory that overcomes 
these practical problems. This must be done without compromising a rigorous treatment 
based firmly on the principles of non-equilibrium quantum statistical physics. Only a theory 
of this nature can address the failures of conventional inflaton dynamics. In this paper such 
a theory is developed. The main result is that, for a minimally coupled scalar field in a de 
Sitter phase, the quantum dynamics of the local order parameter <fr can be described by the 
relatively simple stochastic quantum mechanical Hamiltonian 

H(t) = ^m + e mt V{<P) - ^e" m V'imit), (1-5) 

where p = e 3Ht (p is the canonical momentum and p and <j) obey the usual quantum mechanical 
commutator. F c (t) is a zero mean gaussian colored noise of unit amplitude with a correlation 
time of the order H . This result is valid for a general inflaton potential. The origin of the 
noise is the backreaction of quantum fluctuations with wavelengths shorter than the coarse- 
graining scale. The noise is of a multiplicative nature because its origin is the mode-mode 
coupling induced by the self-interaction of the inflaton. For a free field this theory predicts no 
noise term. This should be expected from arguments made here and previously |12|,|T6|,[13|,|15[ . 
The major simplification is made by ignoring information about spatial correlations between 
the order parameters of different regions. This allows a description based on a single degree 
of freedom. Further significant simplification is obtained by invoking the standard slow 
roll assumptions. This makes it possible to show that the potential renormalization and 
non-local dissipation terms are negligible. 

The approximate semiclassical limit of the quantum open system (1.5) is 

• V'(d>) H 1 / 2 

In deriving this equation the slow roll assumptions we used to neglect the inertial term and 
approximate the colored noise by a white noise F w . This equation is no more complicated 
than the widely used Starobinsky equation (1.4). However it fundamentally differs from it 
because the noise term vanishes for a free theory. Equation (1.6) is a superior alternative to 
the Starobinsky equation because it is the semiclassical limit of the quantum open system 
(1.5). This allows the validity of the semiclassical limit to be dynamically derived by inves- 
tigating the quantum decoherence generated by the noise in (1.5) (On the other hand, the 
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'classicality' of the Starobinsky equation is simply postulated and is not a natural outcome 
of any quantum open system). The other crucial advantage is that, as we will see, equation 
(1.6) leads to a dramatic easing of the fine tuning constraints. 



II. STOCHASTIC INFLATION AND QUANTUM BROWNIAN MOTION 



In this section we will see how the dynamics of a coarse-grained scalar field in an expand- 
ing universe, can be described in terms of time dependent quantum Brownian motion with 
non-linear system-bath coupling. We consider a minimally coupled scalar field evolving in 
a spatially flat background with the metric 



The action has the form 



ds 2 = dt 2 - a 2 (t)d 2 *. 



S = f ds ( d 3 x £(x, s) 

Jti Jn(s) 



(2.1) 



(2.2) 



where 



£(x, s) = a 3 (s) 



(y$(x, s))' 



-\/($(x, 



(2.3) 



We assume is our initial time and a(t«) = 1. 

As is usual in inflation, we split $ into a homogeneous piece and an inhomogeneous piece 

$(x, S ) = 0( S )+V(x,s), (2.4) 

and interpret as the field $ coarse-grained over a volume Vl(s). We therefore write 



1 f 

^) = H7T $(x,s)d 3 x 



which will be true as long as 



n(s) 



rf 3 x^(x, s) = / <i 3 x ip(x, s) = 0. 



n(s) 



Assuming <p ^> ip, the potential term in (2.3) becomes 

\/($)^\/(0)+^'(0) + ^\/"(0). 
Substituting (2.7) into (2.3) and using (2.6) we find 

S = J'ds n(s)a 3 (s) (^0 2 (s) - W)) -\j t tds a 3 (s)V"((f>) U ) d 3 x ^ 2 (x, s) 



(2.5) 



(2.6) 



(2.7) 



+ 



f ds f 

Jti Jtt(s) 



d 3 x a 3 (s) 



1 • (V0(x, s 



(2.8) 
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The coarse-graining breaks the scalar field up into cells of volume f2, each labeled by a 
coarse-grained position variable. We have dropped all reference to the position label on 
since the problem is spatially invariant. By dealing with only one degree of freedom we are 
ignoring information about spatial correlations between different cells. This is the reason 
why the middle term in the expansion (2.7) plays no role. Later in this section we justify 
why we can drop the spatial gradient term of the system sector. 
The environment field ip can be written as 



E^k cosk ■ x + g k sink ■ x] 



(2.9) 



in which case the action (2.8) becomes 



S= / dsQ(s)a 3 {s) 



1 • 



ns)-V{4>) 



EE/, a\s)V"{<j>{s))q^ 



a k 



<r k 



«> 2 " ^ 



12.10) 



where k = |k|. The action (2.10) is now in the form of a time-dependent one dimensional 
system, non-linearly coupled to an environment of time dependent harmonic oscillators. This 



action is a generalization of the well known quantum Brownian motion problem [21,20,40,41 



In the next section we will show how the influence functional method can be used to describe 
the averaged effect of the environment on the quantum dynamics of the system. 



A. Low Frequency Cutoff 

We will impose a long wavelength cutoff in the field ifi, such that the longest wavelength 
equals the diameter of the coarse-graining volume. This is necessary since we are interpreting 
cf) as the field $ coarse-grained over the volume Q(s). This means that all frequencies less 
than the cutoff are considered part of the system sector. Since there is no mode-mode 
coupling for free fields, we cannot expect ip to couple to <fi in this case. This is clearly true 
for the action (2.10). Writing the volume element as 

n « = *& (2 - n) 

we can express the low frequency cutoff as 

k min = eA(t) (2.12) 

where e is an arbitrary dimensionless parameter that scales the coarse-graining volume. An 
obvious question is how do we choose k m i n ? 

At a purely technical level we want to choose Q large enough so that the spatial gradient 
term in (2.3) can be neglected for the system sector. The appropriate volume can be deduced 
from the classical equation of motion for the environment modes. Writing = ^/a(s) we 
find this equation is 
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d 2 v Q k + Q k (k 2 -(d 2 v a)/a) =0 



(2.13) 



where f] = J dt/a(t) is conformal time. We can see that the environment spectrum is split, 
with an unstable low frequency sector and a stable effectively flat space sector. Choosing 

A(t) = yj(d*a)/a, e«l (2.14) 

amounts to including all those modes where the spatial gradient is significant in the envi- 
ronment sector. 

Detailed results in this paper will be limited to de-Sitter inflation where 
a( s ) = e Hs , T] = -e~ Hs /H and U = 0. In this case (2.14) and (2.12) lead to 

k min = eHa(t), e < 1. (2.15) 

This cutoff satisfies our technical requirement. We also have the physical requirement that 
information about large scale structure is contained in the system sector. This requires e 
to be small but finite. Taking e — ■> would correspond to an over coarse-grained system 
that would contain no information of observable scales. The scale of the observable universe 
today is expected to cross the Hubble radius at least 7 e- foldings after inflation begins [|TJ. 
This means that todays horizon scale maps to k ~ 10 3 H. We can see from (2.15) that after 
a few e-foldings, scales corresponding to our observable universe are shifted into the system 
sector. Of course, (2.15) is the same as Starobinsky's choice |5|,|| for the low frequency 
cutoff. 



III. INFLUENCE FUNCTIONAL 

Consider the quantum system described by the action 

S[<t>,q\ = S[<i>] + S e [q\ + SU0,q]. (3.1) 

We will take as our system variable and q to be our environmental variables. Typically the 
environment has infinite degrees of freedom, denoted here by a bold type. We will briefly 
review here the Feynman- Vernon influence functional method for deriving the evolution 
operator [ETj . The method provides an easy way to obtain a functional representation for 



the evolution operator J r of the reduced density matrix p r . 
We are interested in the reduced density matrix defined as 



+oo +oo 



p r (<i>,<f>')= J dqj dqV&q^'.q'Wq-q'). (3.2) 

— oo — oo 

which is propagated in time by the evolution operator J r 

+oo +oo 

p r ((f),<j)',t) = / d<f)i / tfy'i J r {<f>,<l/,t | ()r{<l>i,<f>i,ti ). (3.3) 



By using the functional representation of the full density matrix evolution operator given 
in (3.2), we can also represent J r in path integral form. In general, the expression is very 
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complicated since the evolution operator J r depends on the initial state. If we assume that 
at a given initial time ij, the system and the environment are uncorrelated 



P(ti) = Ps(ti) x p e (ti 



(3.4) 



then the evolution operator for the reduced density matrix does not depend on the initial 
state of the system and can be written as 



Jr(<j>f,<Pf,t I <t>h<t>'vU) 



4> f *i 



' exp {i{S[<f)} - Sm}} m<P'}- (3-5) 



The factor !F[(j), <p'], called the 'influence functional', is defined as 

JF[0, (/)'} = exp {iS IF [4>, <p']} 

= Tr(u[</Hfy]pe{ti)fr[<li 



(3.6) 



where U is the quantum propagator for the action S e [q\ + Si nt [4>(s), q] with 0(s) treated as 
a time dependent classical forcing term. We have found this form to be very convenient for 
deriving the influence functional [ 40|, 41]] . S/^[0, 0'] is the influence action and the effective 
action for the open quantum system is defined as S e ff[4>, 0'] = S[4>] — S[4>'} + Sjf[4>, <fi'}. 

If the interaction term is zero, then it is obvious from its definition, that the influence 
functional is equal to unity and the influence action is zero. In general, the influence func- 
tional is a highly non-local object. Not only does it depend on the time history, but -and 
this is the more important property- it also irreducibly mixes the two sets of histories in the 
path integral of (3.5). In those cases where the initial decoupling condition (3.4) is satisfied, 
the influence functional depends only on the initial state of the environment. The influence 
functional method can be extended to more general conditions, such as thermal equilibrium 
between the system and the environment [f^] , or correlated initial states |21 . 



The influence action for the open system (2.10), in a de Sitter phase, is derived in 
appendix A. It has the form 



S 



IF 



f ds A(s)f(s)+ f ds f ds' A(s)£(s>(s,s') 
Jo Jo Jo 



+ i ds ds' A(s)A(s>(s,s') 
Jo Jo 



where the sum and difference coordinates are 



E(s) = Uv'\<P{s)) + VW(*))), A( S ) = V'Ma)) - V"(<f>'(s)). 



(3.7) 



(3.8) 



The functions v and /i are known respectively as the noise and dissipation kernels. This 
is because, as will be shown in section IV, they correspond to colored noise and non-local 
dissipation. The function f(s) represents a renormalization of the potential. These functions 
are derived exactly in appendix A. 

The approximate form of the noise kernel derived in (A44) is 
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v(s,s')~ cos(4esinh(ifV 2 )) +0(e 2 ), e < 1 and 5 < t c (3.9) 

where 

S = s - S '^ a = (s + s')/2. (3.10) 

As discussed in appendix A, r c is the correlation time of the noise kernel which is defined 
as the time when the argument of the cosine equals 7r/4. The correlation timescale then 
becomes 

t c = 2#- 1 ln(7r/8e). (3.11) 

The exact expression for the dissipation kernel is derived in (A46). We do not quote it 
here since it is a complicated expression that will play little further role. Inspection of 
(A41) and (A46), the exact noise and dissipation kernels, show that they quickly becomes 
highly oscillatory when 5 is greater than the correlation time r c . This means the noise 
and dissipation kernels are effectively cutoff beyond the correlation time. This oscillatory 
behaviour is an artifact of the discrete low-frequency cutoff used for the bath. From (A36) 
the potential renormalization term is 

/(^)^^e 3 ^(l/2 + ln(l/e)), e « 1. (3.12) 

One of the most attractive features of this theory is that the approximate noise and dissi- 
pation kernels (3.9) and (A48), along with (3.12) have no algebraic dependence on e. This 
means that physical results will only very weakly depend on e. 



A. Local Approximation to the Influence Functional 

We can neglect the potential renormalization term in the influence action (3.7) when 

m>>f^' , (0)ln(l/e). (3.13) 

This holds when the potential satisfies the standard 'slow roll' condition Jl| 

\V"((f))\ < 9H 2 ~ 24irV((j))/m 2 pl . (3.14) 

It also holds more generally as long as H -C m p i. 

The slow roll conditions of inflation are necessary to ensure that an approximate de-Sitter 
inflation phase can proceed. In the slow roll limit the dynamical timescale of <fi is much longer 
than the Hubble timescale, so we expect to be able to make a local approximation to the 
dynamics. We start by observing that the integrals over s' in the influence action (3.7) are 
strongly suppressed for 5 > r c . Therefore the lower limit of the integral for s' can be cutoff 
at s — r c . We can now invoke the slow roll condition. On the time scale H^ 1 we expect 
V"((p) to change little, so we can pull out A(s') and S(s') from the integral over s' in the 
influence action. In this case the influence action becomes 
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S IF [<j>,4/\ ~ J ds A(s)£(s) J ds' n(s,s')+i J ds A 2 (s) J ds u(s,s'). (3.15) 

S-T c S — T c 

For the noise integral we now find 

7 v ' ; 647T 6 J s -r c 



# 3 



1927T 6 

where we used 



(3.16) 



cosrr ~ I. (3-17) 
For the dissipation kernel we use the result (A52) derived in appendix A which is 

/ ds'fi(s,s') ~ — — ln(7r/8e). (3.18) 

S-T c 

The influence action (3.15) now becomes 
S 1F \m - J ds e^{v"\<P(s)) - V"\<t>'{s))) + ^jds e^A%s). (3.19) 



We see that the local approximation to the dissipation kernel has reduced to an effective 
potential term in the influence action. This term will be much smaller than the system 
potential when 

V" 2 U) 

V{<t>) » ln(7r/8e)-gAp. (3.20) 

This will hold whenever the slow roll condition (3.14) is satisfied. As an example we can 
substitute V = X(f> 4 into (3.20) which leads to the condition 

A<< l8«' (3 ' 21) 

This is clearly satisfied if A <C 1 which is the standard weak coupling condition. Neglecting 
the potential term in (3.19), the effective action becomes 



S eff [<f>, 0'] = S[4] - SW} + ^ jf* ds e 6 ^A 2 ( s ) (3.22) 

with S[4>] our scaled system action given by (A33) for the case of a de Sitter expansion 
phase. 

In general one must be careful in neglecting non-local dissipation in quantum open sys- 
tems. Dissipation and noise, connected via a fluctuation-dissipation relation, are twin phys- 
ical processes that reflect the relation between energy flowing into and out of the system. 
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This relation is especially important when one is interested in possible stationary states 
where a balance is eventually reached. In inflation we have a very flat potential well away 
from its minimum, and we are only interested in the dynamics over some relatively small 
finite time. As we showed above it is these very special slow roll conditions of inflation that 
allowed us to neglect dissipation. This will certainly not be possible during the reheating 
phase but it should be a good first order approximation to the dynamics in the early slow 
roll phase. 



IV. STOCHASTIC INTERPRETATION OF THE INFLUENCE FUNCTIONAL 



In this section we will show that a stochastic forcing term will generate the imaginary 
part of the influence action (3.7). Consider the action 



S[<j>(sU(s)] = fds (L(0, j>, s) + V"(0)£(s)) 

J ti 



(4.1) 



where £(s) is a zero-mean gaussian stochastic force. This system generates the influence 
functional 



F[$,4T\ = (exp ij*ds A(a)£(s) 



(4.2) 



where the average is understood as a functional integral over £(s), weighted by a normalized 
gaussian probability density functional P[£(s)]. The averaging can be performed to give [H 



t s 

^[0,0'] =exp|-y" ds J ds A(s)A(s>(s,s')j 



(4.3) 



tj 



where u(s, s') is the second cumulant of the force £. Clearly then, before averaging over the 
noise, we can write our effective action derived from (3.7) as 



S eff [<!>, </>', £] = S[<j>] - S[<j>'} + J ds A(s) J ds' X(s')pl( S , s') - f(s) + £(s 



(4.4) 



This result can be used to show that that the dynamics of the effective action (3.22), is 
equivalent to that generated by the quantum mechanical stochastic Hamiltonian 



H{t) 



„2 rr3/2 



2e 3Ht 



+ e 6m V{4>) - 



96tt s 



V"(<j>)F w (t), 



(4.5) 



where p = e 3Ht <f> is the canonical momentum and p and <fi obey the usual quantum mechanical 
commutator. F w (t) is a zero mean gaussian white noise with the correlation function 



(F w (t)F w (t')) = 5(t-t f ). 
The white noise Hamiltonian (4.5) is an approximation to 



(4.6) 
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Hit) = ^btt + ^ Ht V^) - — 3 e 3m V"(cf>)F c (t), (4.7) 

where F c (t) is a zero mean gaussian colored noise of unit amplitude with the correlation 
function 

(F c (t)F c (t')) = cos(4e sinh(if 5/2)) . (4.8) 

This follows because the imaginary part of (3.22) is an approximation of the imaginary part 
of (3.7). Equations (4.5) and (4.7) are the main results of this paper. In the semiclassical 
limit the system timescale is clearly much greater than the Hubble time due to the slow 
rolling nature of inflation. In this case a white noise approximation should be reasonable. In 
the quantum regime it is not so clear that we can approximate (4.7) with (4.5). Numerical 
simulations will be required to answer this. 



A. Semiclassical Limit 

The semiclassical equation of motion for the general effective action (4.4) can be derived 
using 



<,(s e// [0,0',£] 



0, (4.9) 



where = — 0'. With the system action (A33) it can be shown that this equation 
becomes 

+ 3°0 + V (</>)) a 3 (t) + V"'(s)f(s) - V m (4>) f iM(t, s)V"(<f>(s))ds = V""(0)£(t), (4.10) 

where (£(£)£(£')) = v{t,t'). Under special circumstances \i tends to the derivative of a delta 
function which generates local dissipation, hence \x is referred to as the dissipation kernel. 
More generally we see that in the semiclassical limit \x generates non-local dissipation. For 
the influence action (3.7), the functions / and /x are unimportant, and the semiclassical 
equation (4.10) becomes 

+ 3H( p + v'tf) = ^ V "'^) Fc{ t) (4.11) 

where F c is the zero mean gaussian colored noise with correlator (4.8). This equation is 
the semiclassical limit of the quantum open system whose dynamics is described by (4.7). 
Invoking the slow roll approximation (neglecting the 0) this equation becomes 

<i ' +Y w- = ^ v '"^ FM - (412) 

The neglect of the inertial term in (4.11) is standard practice in inflation when describing 
the dynamics of the mean field <f> m . However it should also be a good approximation for 
describing the fluctuations 50 about the mean field (0 m is the solution of (1.3)). Because 
these fluctuations are small we can linearize (4.11) about m to obtain 
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8cj> + 3H8<j> + V"{<j> m )5(f> = ^-V'"{<t> m )F c (t). (4.13) 

The left hand side of this equation describes a damped simple harmonic oscillator. The 
standard over-damping condition is the same as the slow roll condition (3.14). Therefore 
we can neglect the inertial term in (4.13) and conclude that (4.12) should also be a good 
description for the fluctuations as well as the mean field. 

The slow roll conditions allowed us to derive the local approximation (3.22) to the influ- 
ence functional (3.7). In this approximation (4.12) becomes 

where F w is a zero mean gaussian white noise with correlator (4.6). This equation is the 
semiclassical limit of the result (4.5). It is a superior alternative to the Starobinsky equation 
(1.4) for reasons that are addressed in the discussion and conclusion section. Since this 
white noise equation is an approximation to the colored noise equation (4.12), it should be 
interpreted in the Stratonovich sense. 

We have motivated the relatively simple description (4.14) by invoking the slow roll 
conditions and assuming that a semiclassical description is valid. Ultimately the justification 
of equation (4.14) will require a detailed study of the quantum dynamics described by the 
Hamiltonian (4.7). The complicated correlation function (4.8) would make this a difficult 
task. However it must be emphasized that the complex oscillatory nature of this correlation 
function is purely an artifact of the discrete low frequency cutoff used for the bath. The only 
essential 'invariant' information in (4.8) is the correlation time (3.11). Therefore it would be 
reasonable to investigate the effects of colored noise by replacing the correlator (4.8) with an 
exponential correlation function with the correlation time (3.11). This correlator generates 
the simplest type of colored noise dynamics [E3 and would make the quantum dynamics 
described by (4.7) a tractable problem. 



V. AMPLITUDE OF DENSITY FLUCTUATIONS AND THE FINE TUNING 

PROBLEM 

In this section we will use equation (4.14) to estimate the amplitude of primordial density 
fluctuations generated by scalar field fluctuations in inflation. We will consider the potential 
V = A0 4 , for which the solution to the slow roll equation (1.3) is 



0(t) = o exp \-\ — rripit . (5.1) 




With this solution it can be shown that the two slow roll conditions of inflation |TJ break 
down for ~ m p i. We therefore require 0o ^ m pi- The number of e- folds in the inflationary 
era is given by |]J 

8tt ph V{4>) (f>l 
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A solution to the horizon problem requires iV > 60. This means we require 4>o to be at least 
several times the Planck mass. 

At this point it is useful to compare the dynamical timescale and the Hubble timescale. 
The dynamical timescale Td is found from the solution (5.1) to be 




3tt i 

T d = \l — m pl ■ I 5 ' 3 ) 



Using 



H\<P) = f^V(<P) (5.4) 
3m 



P l 



we find that 

Td 



H 



27T0 2 /mJ. (5.5) 



When <p ^> rripi the system timescale is much greater than the Hubble timescale. This 
provides a explicit example of the condition we assumed to make the local approximation 
to the influence functional. 

So far we have only constrained 0o to have a lower bound. In chaotic inflation we expect 
that initially V ~ mi. This means that that we can only make 0o ^ m pi if we a ls° 
correspondingly make A -C 1. In inflation it is the value of the density perturbations which 
serve to fix a value for A, which by the potential constraint also fixes fa. 

The density contrast Sp/p, that is generated by scalar field fluctuations 5<p, is approxi- 
mately given by |L]|| 

6 -P ~ Sr 6<J>. (5.6) 
P <p 

Assuming 3> 8<f), we can linearize equation (4.14) to obtain 

V"(d>) H 1/2 

We can further simplify this equation by putting <p ~ m p i which corresponds to the end of 
inflation. From equation (5.4) the Hubble constant now becomes 

H 2 = yAmJ. (5.8) 

Equation (5.7) is now an Ornstein-Uhlenbeck process for which the stationary distribution 
gives 

-.x 2X2m ll , X 

) - np*- (^) 



In this case (5.6) becomes 
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^ ~ A. (5.10) 

Observational constraints then require A ~ 10~ 6 . Performing an exactly analogous calcula- 
tion for the Starobinsky equation (1.4) gives 

^~v/A (5.11) 
P 

which leads to the conventional result A ~ 10~ 12 . 

Clearly the theory developed here leads to a dramatic easing of the fine-tuning required 
in this model. The result (5.10) was obtained previously by Calzetta and Hu [15j using 



an entirely different approach based on coarse-graining graviton degrees of freedom. They 
pointed out that this new value of A is consistent with the inflaton taking part in nonabelian 
gauge interactions with a coupling constant of 10~ 2 . 



VI. DISCUSSION AND CONCLUSION 

In this paper we have argued that the conventional theory of inflaton dynamics IJ, along 
with Starobinsky's theory of stochastic inflation were problematic because they could 
not explain the quantum-to-classical transition and they predicted an overproduction of 
primordial density fluctuations. It was concluded that a theory which could address these 
issues would involve a conceptual shift that comes with being based on the principles of 
non-equilibrium quantum statistical physics. There have been initial attempts to do this 
[16|,19] using the system environment split originally proposed by Starobinsky ||. However, 



as pointed out in the introduction, these models are very complex and this has so far 
prevented them from being used in a detailed study of inflaton dynamics. 

In this paper a new simplified theory of inflaton dynamics was developed. This ap- 
proach is free of the problems associated with the conventional theory. It is based on the 
Hamiltonian (4.7) which describes the quantum dynamics of the inflaton which has been 
coarse-grained over a constant proper volume of a de Sitter inflationary phase. The origin 
of the noise source is the backreaction of quantum fluctuations with wavelengths shorter 
than the coarse-graining scale. The noise is of a multiplicative nature because its origin is 
the mode-mode coupling induced by the self-interaction of the inflaton. For a free field this 
theory predicts no noise term as we should expect [12,16,13,15||. Apart from addressing the 



problems of the conventional theory, the other essential feature is the relative simplicity of 
the results (4.7) and (4.14). The major simplification was made by ignoring information 
about spatial correlations between the order parameters of different regions. This allows a 
description based on a single degree of freedom. Further significant simplification was ob- 
tained by invoking the standard slow roll assumptions. With this we were able to argue that 
the potential renormalization and non-local dissipation terms were negligible, and that the 
colored noise could be approximated by a white noise in the semiclassical limit. In this case 
the theory is independent of the arbitrary parameter e. The correlation time does depend 
on e but only logarithmically. The robustness of this theory to the coarse-graining volume, 
which is parameterized by e, must be considered a major success. Another very important 
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feature for applications to inflation is that the theory describes a minimally coupled scalar 
field with an arbitrary potential. 

The slow roll semiclassical equation of motion (4.14), is a superior alternative to the 
Starobinsky equation for a number of reasons. It is the approximate semiclassical limit of 
the quantum open system (4.7). This allows the validity of the semiclassical limit to be 
dynamically derived. The 'classicality' of the Starobinsky equation is simply postulated and 
is not a natural outcome of any quantum open system. A quantum mechanical descrip- 
tion, for which the Starobinsky equation is the semiclassical limit, has been attempted using 
methods developed some time ago to describe the quantum dynamics of classically dissipa- 
tive systems ||10|| . However this amounts to simply 'quantizing' an unjustified semiclassical 
equation of motion rather than deriving the quantum theory from fundamental physics. In 
any case, such a procedure is highly suspect because it misconstrues the red-shift term 3H<j) 
in the Starobinsky equation as damping. This repeats an old error in quantum mechanics of 
confusing a time dependent mass with a real dissipation |TTJ. The other major advantage is 



that equation (4.14) leads to a dramatic easing of the fine tuning required in the theory. We 
showed this for the \<f) A potential. The fine tuning problem is systematic to all potentials 
of interest in inflation. Similarly, we should expect that equation (4.14) will lead to a large 
easing of the fine tuning required for these potentials. Finally, the dramatic easing of the 
fine tuning constraints will make equation (4.14) much simpler to numerically simulate than 
the Starobinsky equation. Yi and Vishniac (1993) 0] numerically simulated the Starobin- 
sky equation for the A0 4 potential. They used A = 10~ 6 in their simulations because the 
supposedly correct value of 10~ 12 was numerically intractable. The former value of A is a 
natural outcome of the theory presented here. 

As for future work, the quantum decoherence generated by the noise term in (4.7) will 
allow a dynamical investigation of entropy generation and the quantum-to-classical transi- 
tion. Decoherence studies have been previously attempted fi~6yi9f . However, because these 
models were very complicated the authors were restricted to numerically calculating the 
diffusion coefficient in the quantum master equation. Unfortunately this tells us nothing 
about whether or not enough decoherence will occur in the model to justify a semiclassical 
description. The critical dynamical aspect to this problem is the competition between the 
coherence generating non-linear potential and inflation expansion effects, and the decoher- 
ing stochastic term. This can only be investigated with a full solution of the non-unitary 
quantum dynamics. If the noise term is strong enough we would hope that decoherence will 
justify a semiclassical description. However, we must also remember that the noise can not 
be so large as to destroy the essentially deterministic slow roll dynamics of the local order 
parameter 0. The inflaton dynamics described by equation (4.7) is the first viable theory to 
permit a detailed study of the quantum-to-classical of the local order parameter <fi. This is 
the critical process that leads to the generation of classical density fluctuations. 

Assuming that the semiclassical description (4.14) is valid, we can use this equation to 
deduce the amplitude and spectrum of the generated primordial density fluctuations, along 
with the implications this has for coupling constants of various inflaton potentials. Also 
of great interest is the generation of any relevant non-gaussian features in the resulting 
probability distribution of <fi 0, and the implications of equation (4.14) for the very large 
scale structure of the universe ||. Till now, these issues have been addressed with the 
Starobinsky equation as their theoretical foundation. An interesting technical issue is what 
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effect relaxing the white noise approximation would have on these results. Investigations 
into the implications of the theory presented here are in progress. 
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APPENDIX A: CALCULATING THE INFLUENCE FUNCTIONAL 



The exact influence functional for the action of the type (2.10) has been previously found 
to be 141 



(Al) 



where a, (3 are Bogolubov coefficients which must satisfy the constraint 

\a\ 2 - \(3\ 2 = I. 



(A2) 



The influence functional (Al) differs from that in |^T|] because we have here taken into 
account the double set of modes in the action (2.10). The Bogolubov coefficients are solutions 
of the equations 



where 



9W 



a = —ig(3 — iha 
(3 = ih(3 + iga 



2 I k 



a 3 (s) 



(A3) 
(A4) 



(A5) 



* (a) = 2bw + -nb— J' (A6) 

and 

oo 2 (s) = k 2 /a 2 + V"(<f)(s)). (A7) 

We must have a = 1 and (3 — at U as our initial conditions. 
Defining 

U(M<) = f *H ) • (A8) 



Pl<f>t,„] C\4, 



t . t , 



the solution of (A3-4) can be written as 
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where 



We will write 



where 



U(t, U) = T exp (-i £ ds u(s)) (A9) 



(A10) 



u(s) = u (s) + ui(s) (All) 



.w-M^M. (A12) 



In this case we have 



U(*,t i )=U (t,t i )Ui(t,t i ) (A13) 
where U (t,tj) is U(t, tj) evaluated at V"(0) = 0. The equation of motion for Ui is 

Ui(s,ti) = -2(Uo 1 (s,t t )ui(s)U (s,t,))Ui(s,t i ) (A14) 
where U -1 denotes the inverse matrix. This inverse is 



u-\t,t l ) = ( ° 



because of the constraint (A2). From (A13) and (A8) we have 

a[<jH,u\ = <xo(t,t i )a 1 [<iH,t i \+ /3 *(Mi)/Wt,*J (A16) 

=/3 o (t,t l )ai[0 M J+a*(t,t l )/3 1 [0 M J. (A17) 
Substituting this into the influence functional (Al), and using the constraint (A2), we find 

n<f>,<t>'\ = iiW^km - PiWmy 1 (Ais) 

k 

where we have dropped the k subscript on the Bogolubov coefficients. 
To second order in V" the solution to (A14) is 

rt 



V 1 (t,t l ) = l-i ( ds U 1 (s,t l )u 1 (s)U (s,t J ) 

- f ds f ds' U 1 ( S ,t J )ui(s)Uo(s,t l )U 1 ( S , ,t,)u 1 ( S / )U (s , ,t i ). (A19) 



Using this solution and (A18), the influence action becomes 
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S IF [<f>,<ff\ = Y,{-^f t ds a 3 (s)A(s)X k (s)X* k (s) 



+ W 2 Ju dS L dS ' a3 ( s ) a3 ( s ')A( S )A( S ')K( S )X fe 2 ( S ') + X 2 (s)Xf (s>) 
~ht dS \l d * (s)Xl(s') -X 2 k (s)X* k 2 (*')]} (A20) 



where 



and 



A(s) = V"M8)) - V"{<t>>{s)), E(s) = - [V"{<t>{s)) + V"{<t>'{s)) 



X k (s) = a (s, U) + Pq(s, U). 

Using (A3-4) we can show that (A22) obeys the classical equation of motion 

a • k 2 
X + 3-X + — X = 

a cr 

subject to the initial conditions 

x(u) = i, x(u) = -ik. 

These initial conditions ensure a(ti) = 1 and /3(ti) = 0. 
Taking the continuum limit we have 

£ - 2L [ k 2 dk. 

k 



4tt 2 Jk. 



(A21) 



(A22) 



(A23) 



(A24) 



(A25) 



where we adopt the notation that Q is a symmetric function of s and s', i.e. Q 2 = Q(s)Q(s'). 
In this case the influence action (A20) becomes 



S IF [(f),(f)'}= - f dsA(s)f(s) + f ds f S ds' A(s)S(s')/i(s,s') 

+ i f ds f S ds' A(s) A(s>(s, s') (A26) 



where 



f(s) = a 3 (s) 



8^ J k 



dk kX k (s)X* k (s), 



(A27) 



u{s,s') = a 3 {s)a 3 (s') 



n 



and 



32tt 2 



Q roc r - 

H(s,s') = -ia 3 (s)a 3 (s')-— dk [xf(s)X 2 k (s')-X 2 k (s)Xf(s')}. (A29) 



dk 



Xf(s)X 2 k (s')+X 2 k (s)Xf(s') 



(A28) 
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From the action (2.8) we see that it is ip 2 that couples to the system. The truncation 
of the solution (A19) amounts to assuming that ip 2 generates a gaussian noise source on 
the system. This is proved in section 4. The gaussian property ensures that only the first 
two cumulants of the stochastic process generated by ip 2 appears in the influence functional. 



The influence functional can be expressed in the alternative form [IT 



/i(3,/)~^ 2 (^ 2 ( s '))-^ 2 W 2 ^)) 

u(s, s') ~ ty 2 (a)VV)> + <V> 2 (*V(s)> - 2<^( S )><VV)> 



(A30) 



which makes this clearer. For a gaussian process all higher order cumulants vanish because 
all higher order moments can be expressed in terms of first and second moments. Clearly 
then, the truncation of the solution (A19) is an equivalent approximation to writing all 
higher order moments of if) 2 in terms of first and second moments. 
Our quantum field theory has the commutator 



*(x),P*(y) =iM(x-y) 



(A31) 



where P$(y) is the canonical field momentum derived from the action (2.3). Using this 
commutator we can show that our coarse-grained field (2.5), and a similarly defined coarse- 
grained canonical momentum, will obey the usual single particle quantum mechanical com- 
mutator only after we perform the scaling 

ft -> Q(s)h (A32) 

on the Planck constant. This is a natural consequence of reducing a quantum field theoretic 



problem to a quantum mechanical problem and has been discussed previously JTTJ. We have 
set ft = 1 in this paper, but the scaling (A32) is equivalent to scaling the effective action 
(4.4) by Q(s). This gives the new system action 



ds a 3 (s) 



v(<p) 



and it requires that we must scale the functions (A27-A29) as 



(A33) 



(A34) 



Equation (A23) is identical to (2.13) when we scale X by the scale factor and transform 
to conformal time. For a de Sitter phase, the solution of (2.13) leads to 



X k (s) 




(A35) 



Strictly speaking an appropriate linear combination of these complex mode functions are 
required to satisfy the initial conditions (A24). However the time dependent low frequency 
cutoff ensures that any dependence on the initial condition is a transient effect. So the 
description based on the mode function (A35) is accurate on the Hubble timescale. With 
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(A35) we are now in a position to calculate (A27-A29) using (2.15) as the low frequency 
cutoff. 

We will first consider the function f(s). This function is ultraviolet divergent which is not 
surprising since / ~ (4> 2 )- We will adopt the usual non-rigorous renormalization procedure 
here simply by imposing the inflationary ultraviolet cutoff k = Ha (see Brandenberger 
(1984) 0] and Habib [14]|). This cutoff is determined by demanding that the energy density 
of inflaton fluctuations be less than the potential (which is supposed to dominate), and that 
the inflaton fluctuations can be taken to be in the vacuum state. This is discussed by Liddle 
and Lyth JTJ. We must also perform the scaling (A34) on (A27) which eliminates the volume 
factor Q from (A27). We therefore find 



/(«) 



a(s) 



8tt 2< 



Ha 1 

dk\k + — 



eHa 



^3Hs 



krf 

l/2 + ln(l/e)), e<l. 



(A36) 



The small e limit is necessary in order to ignore the spatial gradient term in the system 
sector. This was shown in section 2 where it was also discussed that we require e to be 
small but finite. We are therefore never faced with any problems associated with infrared 
divergences. 

After performing the scaling (A34) on (A28) and substituting (A35) into (A28) we find 
that the noise kernel becomes 




^ 2 
'/ 



1 



k 2 r]7]' k 2r q 2r q' 2 k 4 rfr]' 



/2 



+ 2<L sm2kS, 



where 



7] — 7] 



(A37) 



(A38) 



The noise kernel is clearly ultraviolet finite and therefore requires no renormalization. This 
should not be surprising. From (A30) we know that the kernels are built from the quadratic 
field operators ip 2 rather than the linear operators ip of conventional two-point functions 
(which are well known to be ultraviolet divergent). The integrals in (A37) can be done by 
using the integral identities 



ikx 



dk 



1 



Akx 



„m— 1 



IX 



Akx 



~m—l 



dk 



(A39) 



Akx 



J^— dk = Ci(kx) + iSi(kx) (A40) 
where Si and Ci are the sine and cosine integral functions . We then find that (A37) becomes 
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rr4 

v(s, s') = J^E e6Ha {^H-^d^) + cos(4e smh(H5/2)) (2 + 12e 2 + 8e 2 sinh 2 (iJ5/2)) 
+ sin(4esinh(iJ5/2)) (^8esmh(H5/2) - ^ smh- 1 (H5/2)j 

+ 16e 3 Si(4e smh(H5/2)) (3 sinh(if 5/2) + 2 sinh 3 (iJ5/2))] , (A41) 

where 5d denotes the Dirac delta function which is not to be confused with 8 defined in 
(3.10). This noise kernel does not diverge for f < 1 because the scaling (A34) makes the 
kernel proportional to fi _1 . Using (2.11) with A = Ha we find that 

o 3 fj3 p 3Hs 

= (A42) 

which shows that fi -1 scales the kernel by e 3 . The noise kernel (A41) has a correlation time 
which we will define as the time when the argument of the cosine in (A41) equals 7r/4. We 
then find that the correlation timescale r c is 

r c = 2#- 1 ln(7r/8e). (A43) 

The noise kernel is highly oscillatory for S > r c . It is therefore effectively cutoff beyond the 
correlation time. We can therefore consider the approximate noise kernel 

ttA 

v(s,s') ~ ^^e 6Ha cos(4esinh(iJ5/2)) +0(e 2 ), e < 1 and 5 < r c . (A44) 

Clearly the noise kernel has simplified greatly. However what is most pleasing is that the 
noise kernel has lost its algebraic dependence on e. Because the noise kernel is ultraviolet 
finite we didn't impose the inflationary cutoff. Had we done so the cutoff would have 
contributed terms of order e 3 into the noise kernel. Thus in the small e limit the noise kernel 
is also independent of any ultraviolet cutoff. 

We will now consider the dissipation kernel (A29). The scaling (A34) eliminates the 
volume factor Q from the kernel in (A29) which upon substituting (A35) becomes 



a(s)a(s') r H " 

1 I 



' / _2_ 8 2 n J_ 

v \ k 2 r]r]' k 2 r] 2 rj' 2 k^rfrj' 2 



25 v cos 2k5 v 



krjrj' k 3 i] 2 r]' 2 



(A45) 



All of the terms in the above integral are ultraviolet finite except for the first flat space 
contribution which leads to a logarithmic divergence. For this reason we have used the 
inflationary ultraviolet cutoff. Using the integrals (A39-40) we find that (A45) becomes 



TT 

3Ha r sin(4e smh(#<5/2)) (2e~ 3 + 12c" 1 + Se" 1 sinh 2 (F5/2)) 



+ cos(4esinh(#(V2)) (J sinh" 1 ^/ 2 ) - 8e~ 2 sinh(F5/2) 
- 16Ci(4esinh(iJ5/2))(3sinh(iJ5/2) + 2sinh 3 (iJ5/2)) 
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- sin(4 smh(HS/2)) (l4 + 8 sinh 2 (/J5/2)) 

- cos(4sinh(/75/2)) (J smh- 1 {H5/2) - 8 swh{H6/2) 
+ 16Ci(4 sinh(iJ5/2)) (3 sinh(if5/2) + 2 sinh 3 (/J5/2))] . (A46) 

We are only interested in the dissipation kernel for 5 < r c . For times greater than this the 
dissipation kernel is highly oscillatory and therefore effectively cutoff. In this regime we can 
use the approximations 

cosrr ~ I — x 2 /2, sinx^x — x 3 /3 (A47) 
to get the leading order approximation to the dissipation kernel (A46) which is 



/i(s, s') ~ ^e 3Ha smh~\H5/2) + sinh(if<5/2) + ^ smh 3 (H5/2) 
n z \_32 9 

2 



- Ci(4esinh(//5/2))(sinh(//5/2) + - sinh 3 (i/5/2) 

7 I 

- sin (4 sinh(if 5/2)) (— + - sinh 2 (iJ5/2)) 

- cos(4 sinh(if 5/2)) sinlr^if 5/2) - ^ sinh(if 5/2)) (A48) 
+ Ci(4 sinh(#5/2)) (sinh(#5/2) + - sinh 3 (#5/2)) + 0(e 2 ) , e< 1 and 5 < r c . 

Remarkably we find that in this approximation the dissipation kernel has lost its algebraic 
dependence on e. 

Unfortunately the dissipation kernel can not be approximated to a simple expression 
like the noise kernel. However the slow roll conditions of inflation ensure that we need only 
consider an order of magnitude estimate of the integral 

/ ds'fi(s, s') ~ r c /i(s, s - t c ) (A49) 

JS-T C 

which we obtain by multiplying the correlation time r c defined in (A43) by the dissipation 
kernel evaluated at 5 = r c . From (A48) we can write 

^ S - Tc ) ~ l^L e 3Hs e -3Hr c /2 sinh ^ Hr J 2 ) ~ JL e ™s (A5Q) 

where we used 

sinh(#r c /2) ~ e HTc/2 /2. (A5I) 
Using (A43) and (A50) we find that (A49) becomes 

f ds'n(s,s') ~ i— ln(7r/8e). (A52) 

Js-T r 4n 2 



25 



REFERENCES 



[1] T. Padmanabhan, Structure Formation in the Early Universe, Cambridge University 
Press (1993); P.J.E. Peebles, Principles of Physical Cosmology, Princeton University 
Press, N.J. (1993); A.R. Liddle and D.H. Lyth, Phys. Rep. 231, 1 (1993); A.D. Linde, 
Particle Physics and Inflationary Cosmology, Harwood Academic Publishers GmbH, 
Chur, Switzerland, (1990); E.W. Kolb and M.S. Turner, The Early Universe, Addison 
Wesley (1990). 

[2] A. Guth and S. Y. Pi, Phys. Rev. Lett. 49, 1110 (1982); A. A. Starobinsky, Phys. Lett. 
117B, 175 (1982); S. W. Hawking, Phys. Lett. 115B, 295 (1982); J. M. Bardeen, P. 
J. Steinhardt and M. S. Turner, Phys. Rev. D28, 629 (1983); R. Brandenberger, R. 
Kahn and W. Press, Phys. Rev. D 28, 1809 (1983); V. Mukhanov, H. Feldman and 
R. Brandenberger, Phys. Rep. 215, 203 (1992); R. Brandenberger, Nucl. Phys. B. 245, 
328 (1984); A. H. Guth and S. Y. Pi, Phys. Rev. D32, 1899 (1985). 

[3] See, e.g., K. Lindenberg and B. J. West, The Nonequilibrium Statistical Mechanics of 
Open and Closed Systems (VCH Press, New York, 1990). 

[4] M. Morikawa, Prog. Theor. Phys 93, 685 (1995). 

[5] A. A. Starobinsky, in Field Theory, Quantum Gravity and Strings, ed. H. J. de Vega 

and N. Sanchez (Springer, Berlin 1986). 
[6] S. J. Rey, Nucl. Phys. B284, 706 (1987); J. M. Bardeen and G. J. Bublik, Class. Quan. 

Grav. 4, 473 (1987); M. Morikawa, Phys. Rev. D42, 1027 (1990); H. E. Kandrup, Phys. 

Rev. 39, 2245 (1989) 

[7] I. Yi and E. T. Vishniac, Phys. Rev. D43, 5295 (1993); Phys. Rev. D45, 3441 (1992); 

D.S. Salopek and J.R. Bond, Phys. Rev. D43, 1005 (1991); I. Yi, E.T. Vishniac and S. 

Mineshige, Phys. Rev. D43, 362 (1991); S. Mollerach, S. Matarrese, A. Ortolan, and F. 

Lucchin, Phys. Rev. D44, 1670 (1991); G.V. Chibisov and Yu. V. Shtanov, Int. J. Mod. 

Phys. A13, 2625 (1990); D.S. Salopek, J.R. Bond and J.M. Bardeen, Phys. Rev. D40, 

1753 (1989); S. Matarrese, A. Ortolan and F. Lucchin, Phys. Rev. D40, 290 (1989); A. 

Ortolan, F. Lucchin and S. Matarrese, Phys. Rev. D38, 462 (1988). 
[8] A. Linde, D. Linde and A. Mezhlumian, Phys. Rev. D 49, 1783 (1994); Phys. Lett. 

B307, 25 (1993); Phys. Lett. B345, 203 (1995); Y. Nambu and M. Sasaki, Phys. Lett. 

B 219, 240 (1989); Y. Nambu, Prog. Theor. Phys 81, 1037 (1989). 
[9] A. Hosoya, M. Morikawa and K. Nakayama, Int. J. Mod. Phys. A 4, 2613 (1989); Y. 

Nambu, Phys. Lett. B 276, 11 (1992); M. Mijic, Phys. Rev. D49, 6434 (1994). 
[10] F.R. Graziani, Phys. Rev. D38, 1122 (1988);D38, 1131 (1988);D38, 1802 (1988); D39, 

3630 (1989);O.E. Buryak, Phys. Rev. D 53, 1763 (1996). 
[11] This is discussed by J. R. Ray, Am. J. Phys. 47, 626 (1979). 

[12] B. L. Hu and Y. Zhang, "Coarse-Graining, Scaling, and Inflation" Univ. Maryland 
Preprint 90-186 (1990); B. L. Hu, in Relativity and Gravitation: Classical and Quantum 
Proceedings of SILARG VII, Cocoyoc, Mexico, Dec. 1990. eds. J. C. D' Olivo et al 
(World Scientific, Singapore, 1991). 

[13] M. Morikawa, in Quantum Physics and the Universe, ed. Namiki, K. Maeda, et al 
(Pergamon Press, Tokyo, 1993). Vistas in Astronomy, 37, 87 (1993). 

[14] S. Habib, Phys. Rev. D46, 2408 (1992) 

[15] E. Calzetta and B.L. Hu, Phys. Rev. D 52, 6770 (1995) 



26 



B. L. Hu, J. P. Paz and Y. Zhang, in The Origin of Structure in the Universe, ed. by 

E. Gunzig and P. Nardone (Kluwer, Dordrecht 1993). 

B. L. Hu, in Quantum Physics and the Universe, ed. Namiki, K. Maeda, et al (Pergamon 
Press, Tokyo, 1993). Vistas in Astronomy 37, 391 (1993). 

B. L. Hu, in Proc. Second International Workshop on Thermal Fields and Their Appli- 
cations, eds. H. Ezawa et al (North-Holland, Amsterdam, 1991). 

F. Lombardo and F.D. Mazzitelli, Phys. Rev. D 53, 2001 (1996). 

B. L. Hu, J. P. Paz and Y. Zhang, Phys. Rev. D45, 2843 (1992); Phys. Rev. D47, 1576 
(1993). 

R. Feynman and F. Vernon, Ann. Phys. (NY) 24, 118 (1963); R. Feynman and A. 
Hibbs, Quantum Mechanics and Path Integrals, (McGraw - Hill, New York, 1965); A. 
O. Caldeira and A. J. Leggett, Physica 121A, 587 (1983); H. Grabert, P. Schramm and 

G. L. Ingold, Phys. Rep. 168, 115 (1988). 
V. Hakim and V. Ambegoakar, Phys. Rev. A 32, 423 (1985). 

B. L. Hu, in Quantum Mechanics in Curved Spacetime, ed. J. Audretsch and V. de 
Sabbata (Plenum, London 1990). 

S. Sinha, Ph. D. Thesis, University of Maryland (1991) 

J. P. Paz and S. Sinha, Phys. Rev. D44, 1038 (1991); ibid D45, 2823 (1992). 
B. L. Hu, J. P. Paz and S. Sinha, in Directions in General Relativity Vol. 1: Misner 
Festschrift eds B. L. Hu, M. P. Ryan and C. V. Vishveswara (Cambridge University 
Press, Cambridge, 1993). 

For earlier work on decoherence in quantum cosmology, see, e.g., C. Kiefer, Clas. Q. 
Grav. 4, 1369 (1987); J. J. Halliwell, Phys. Rev. D39, 2912 (1989); T. Padmanabhan, 
ibid. 2924 (1989). 

B. L. Hu, Physica A158, 399 (1989). 

B. L. Hu, in Proceedings of the Third International Workshop on Thermal Field Theory 
and Applications, ed. R. Kobes and G. Kunstatter (World Scientific, Singapore, 1994). 
B. L. Hu and Y. Zhang, in Proceedings of the Third International Workshop on Quantum 
Nonintegrability, eds J. M. Yuan, D. H. Feng and G. M. Zaslavsky (Gordon and Breach, 
Langhorne, 1993); Mod. Phys. Lett. A8, 3575 (1993). 

J.J. Halliwell, Phys. Rev. D48, 2739 (1993); A. Anderson and J. J. Halliwell, Phys. 
Rev. D48, 2753 (1993); J.J. Halliwell, Phys. Rev. D48, 4785 (1993). 
W. H. Zurek, Phys. Rev. D24, 1516 (1981); D26, 1862 (1982); in Frontiers of Nonequi- 
librium Statistical Physics, ed. G. T. Moore and M. O. Scully (Plenum, N. Y., 1986); 
Physics Today 44, 36 (1991); E. Joos and H. D. Zeh, Z. Phys. B59, 223 (1985); A. 
O. Caldeira and A. J. Leggett, Phys. Rev. A31, 1059 (1985); W. G. Unruh and W. H. 
Zurek, Phys. Rev. D40, 1071 (1989); B. L. Hu, J. P. Paz and Y. Zhang, Phys. Rev. 
D45, 2843 (1992); D47, 1576 (1993); W. H. Zurek, J. P. Paz and S. Habib, Phys. Rev. 
Lett. 47, 1187 (1993); J. P. Paz, S. Habib and W. H. Zurek, Phys. Rev. D47, 488 
(1993); J. P. Paz and W. H. Zurek, Phys. Rev. D48, 2728 (1993). 
[33] R. B. Griffiths, J. Stat. Phys. 36, 219 (1984); R. Omnes, J. Stat Phys. 53, 893, 933, 

957 (1988); Ann. Phys. (NY) 201, 354 (1990); Rev. Mod. Phys. 64, 339 (1992). 
[34] M. Gell-Mann and J. B. Hartle, in Complexity, Entropy and the Physics of Information, 
ed. by W. H. Zurek (Addison- Wesley, Reading, 1990); H. F. Dowker and J. J. Halliwell, 
Phys. Rev. D46, 1580 (1992); A. Albrecht, Phys. Rev. D48, 3768 (1993); J. Twamley, 



27 



Phys. Rev. D48 (1993). 
[35] M. Gell-Mann and J. Hartle, Phys. Rev. D47, 3345 (1993); T. Brun, Phys. Rev. D47, 
3383 (1993). 

[36] E. Calzetta and B. L. Hu, in Directions in General Relativity, Vol II: Brill Festschrift 

eds B. L. Hu and T. A. Jacobson (Cambridge University Press, Cambridge, 1993). 
[37] W. H. Zurek, Prog. Theor. Phys. 89, 281 (1993). 

[38] J. B. Hartle, Les Houches 1992 Lectures (1993). J. B. Hartle, in Directions in General 
Relativity, Vol. 1: Misner Festschrift eds B. L. Hu, M. P. Ryan and C. V. Vishveswara 
(Cambridge Univ., Cambridge, 1993). 

[39] R. Omnes, Rev. Mod. Phys. 64, 339 (1992). 

[40] B. L. Hu and A. Matacz, Phys. Rev. D49, 6612 (1994). 

[41] B. L. Hu and A. Matacz, Phys. Rev. D 51, 1577 (1995). 

[42] L. Grishchuk and Y. V. Sidorov, Phys. Rev. D42, 3414 (1990); A. Matacz, Phys. Rev. 
D49, 788 (1994). 

[43] A.J. McKane, H.C. Luckock and A.J. Bray, Phys. Rev. A 41, 644 (1990); A.J. McKane, 
A.J. Bray and T. Newman, Phys. Rev. A 41, 657 (1990); H.C. Luckock and A.J. 
McKane, Phys. Rev. A 42, 1982 (1990). 



28 



